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The magnetoconductance G in chaotic quantum dots at medium/high magnetic fluxes $ is calcu- 
lated by means of a tight binding Hamiltonian on a square lattice. Chaotic dots are simulated by 
introducing diagonal disorder on surface sites of L x L clusters. It is shown that when the ratio W/L 
is sufficiently large, W being the leads width, G increases steadily showing a maximum at a flux 
<&max tx W. Bulk disordered ballistic cavities (with an amount of impurities proportional to L) does 
not show this effect. On the other hand, for magnetic fluxes larger than that for which the cyclotron 
radius is of the order of L/2, the average magnetoconductance inceases almost linearly with the 
flux with a slope proportional to W 2 , shows a maximum and then decreases stepwise. These results 
closely follow a theory proposed by Beenakker and van Houten to explain the magnetoconductance 
of two point contacts in series. 



I. INTRODUCTION 

Magnetoconductance in chaotic quantum dots have at- 
tracted a great deal of attention in recent years 
Weak localization effects have been thoroughly investi- 
gated searching for differences between chaotic and regu- 
lar cavities [|l]|| . More recently the selfsimilar character 
of magnetoconductance fluctuations in chaotic quantum 
dots has deserved several experimental studies ^Q]- In 
Ref. Q it was reported that in cavities with sufficiently 
soft walls, rather wide leads and a high zero field conduc- 
tance of the order of 40 conductance quanta, fluctuations 
were very weak and the magnetoconductance increased 
steadily in approximately 20% over 50 flux quanta 
Although an increase in the magnetoconductance as a 
function of the magnetic flux in cavities with wide leads 
may not be that surprising, no theoretical analysis of this 
result is yet available. This is so despite of the wealth 
of experimental and theoretical information on magneto- 
conductance in related mesoscopic systems [^[-^0|. 

The purpose of the present work is to investigate the 
magnetoconductance of chaotic quantum dots over a 
wide range of magnetic field. Quantum dots are de- 
scribed by means of a tight-binding Hamiltonian on L x L 
clusters of the square lattice. Non-regular (chaotic) be- 
havior is induced by introducing disorder at surface sites, 
a procedure that has been shown to reproduce all proper- 
ties of quantum chaotic cavities |Tl]] . The most outstand- 
ing conclusions derived from our study are the following. 
For sufficiently open systems, large leads width W or, 
alternatively, high zero field conductance, the magneto- 
conductance increases steadily as a function of the mag- 
netic flux, reaching a maximum at a magnetic flux $ raal 
proportional to the the leads width. This effect, which 
is in agreement with the experimental observations of 



Q, does not show up in regular or bulk disordered cav- 
ities. The average magnetoconductance versus magnetic 
flux curve shows four clearly differentiated regions: i) At 
small fluxes (typically below 1-2 flux quanta) the weak 
localization peak with the typical Lorentzian shape is 
observed, ii) This is followed by a flux range over which 
the magnetoconductance shows a non-universal behav- 
ior which depends on leads configuration, iii) The latter 
lasts until the cyclotron radius becomes of the order of 
L/2. Beyond this point the average magnetoconductance 
increases linearly with the magnetic flux with a slope 
which increases with the leads width (for very small W 
the slope is nearly zero and the magnetoconductance re- 
mains constant in a large flux range) . iv) At large fluxes 
the magnetoconductance decreases stepwise (each step of 
one flux quanta), due to the successive crossings of the 
Fermi energy of the transversal modes that contribute to 
the current. As discussed below these results are com- 
patible with a theory proposed by Beenakker and van 
Houten |l(| to interpret the experimental results for the 
magnetoconductance of two point contacts in series 

The rest of the paper is organized as follows. Section 
II includes a description of our model of chaotic quantum 
dot and of the method we used to compute the current. 
The results are presented in Section III, and discussed in 
terms of the theory of Ref. jl0| in Section III. Section IV 
is devoted to summarize the conclusions of our work. 

II. MODEL AND METHODS 
A. Model 

The quantum dot is described by means of a tight- 
binding Hamiltonian with a single atomic level per lattice 
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site on L x L clusters of the square lattice: 
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where \m, n > represents an atomic orbital on site (m, n). 
Indexes run from 1 to L, and the symbol <> denotes 
that the sum is restricted to nearest-neighbors. Using 
Landau's gauge the hopping integral is given by, 
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where the magnetic flux $ is measured in units of the 
quantum of magnetic flux $o = h/e. The energy uj m . n of 
atomic levels at impurity sites (IS) is randomly chosen 
between —A/2 and A/2, whereas at other sites uj m , n = 0. 
Impurities were taken on all surface sites |ll],[l3|] but those 
coinciding with the leads entrance sites to avoid excessive 
(unphysical) scattering. This model has been proposed 
to simulate cavities with rough boundaries and its prop- 
erties closely follow those characterizing quantum chaotic 
systems pl| , ^3| . Some calculations were also carried out 
on clusters with 2L bulk impurities fll7jl. 



B. Conductance 

The conductance (measured in units of the quantum 
of conductance Go = e 2 /h) was computed by using an 
efficient implementation of Kubo formula. The method 
is described in |l4|] , while applications to mesoscopic sys- 
tems can be found in p5| , ^6[ . For a current propagating 
in the x-direction, the static electrical conductivity is 
given by: 
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where IuyQ(E) is obtained from the advanced and re- 
tarded Green functions: 
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and the velocity (current) operator v x is related to the 
position operator x through the equation of motion Tiv x = 

H,x , H being the Hamiltonian. 

Numerical calculations were carried out connecting 
quantum dots to semiinfinite leads of width W in the 
range 1-L. The hopping integral inside the leads and 
between leads and dot at the contact sites is taken equal 
to that in the quantum dot (ballistic case). Assuming 
the validity of both the one-electron approximation and 



linear response, the exact form of the electric field does 
not change the value of G. An abrupt potential drop at 
one of the two junctions provides the simplest numerical 
implementation of the Kubo formula |l4| since, in this 
case, the velocity operator has finite matrix elements on 
only two adjacent layers and Green functions are just 
needed for this restricted subset of sites. Assuming this 
potential drop to occur at the left contact (Ic) side, the 
velocity operator can be explicitly written as, 
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where (\lc,j > are the atomic orbitals at the left contact 
sites nearest neighbors to the dot. 
Green functions are given by: 
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where T,i^(E) are the selfenergies introduced by the two 
semiinfinite leads ||. The explicit form of the retarded 
selfenergy due to the mode of wavevector k y is: 
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for energies within its band \E—e(k y )\ < 2, where e(k y ) = 
2cos(k y ) is the eigenenergy of the mode k y which is quan- 
tized as k y = (rife tt)/(W+1), being an integer from 1 
to W. The transformation from the normal modes to the 
local tight-binding basis is obtained from the amplitudes 
of the normal modes, < n\k y >= y // 2/(W + l)sin(nfc y ). 
Note that in writing Eq. (0) we assumed that the mag- 
netic field was zero outside the dot ||. This point will 
have some relevance in relation to the interpretation of 
our numerical results in terms of the theory of Ref . (TTJ] . 

Calculations were carried out on clusters of linear size 
L = 47 — 394 (in units of the lattice constant) and at 
a fixed arbitrarily chosen, Fermi energy E = —tt/3. In 
some cases averages over disorder realizations were also 
done. Although most calculations were performed with 
input /output leads of width W connected from site (1,1) 
to site (1, l + I'F), and from (L,L) to (L, L — W), respec- 
tively, other input/output leads configurations were also 
explored. All calculations on disordered systems were 
done at a fixed value of the disorder parameter A=6. 



III. RESULTS 

We first discuss results for a W/L similar to that used 
in the experiments of Ref. [Q . In that work conductance 
measurements were taken on a stadium cavity with a 
lithographic radius of 1.1 /im and leads 0.7 /jm wide, 
which gives a W/L ratio of 0.64. Fig. [j] shows the results 
for the magnetoconductance in cavities of linear sizes in 
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the range L=47-394 and leads of width W = 0.65L. 
The results correspond to a single realization of disor- 
der. The most interesting result is the steady increase 
of the conductance with the magnetic field. The conduc- 
tance reaches a maximum at a magnetic flux which, as 
illustrated in the inset of the Figure, increases linearly 
with the leads width, or the linear size of the system. 
Note that, as the results of Fig. [j] correspond to W oc L, 
they cannot allow to identify which of the two param- 
eters (W and L) control the maximum in G. We have 
checked, however, that it is in fact W the one that mat- 
ters (see also below). The increase in G until the maxi- 
mum is reached can be as high as 30%. It is interesting 
to remark that the increase in the conductance occurs 
with relatively small fluctuations due to the large W/L 
ratio (or degree of opening) of the cavity (see jl6| ) . Al- 
though the experimental data were taken at fields not 
high enough to observe the maximum shown in Fig. [I], 
it can be safely assessed that our results are compatible 
with those of Ref. fl. In particular the zero field con- 
ductance reported in that work was around 38 quanta, 
increasing up to 46 quanta over approximately 50 flux 
quanta. This is rather similar to the magnetoconduc- 
tance curve for L = 197 shown in Fig. |l|. 

In Fig. U we compare the results for the cavity with 
surface disorder with those for a regular cavity and for 
a cavity with 2L bulk impurities. The results indicate 
that the cavity having surface disorder is the only one 
that reproduces the experimental results |0|. In regu- 
lar cavities the conductance does not increase steadily 
due to the large amplitude oscillation discussed in Jl6| . 
This behavior clearly differentiates regular and chaotic 
cavities. At large fields the result for the cavity with sur- 
face disorder coincides with that for the regular cavity. 
This is a consequence of the fact that for sufficiently high 
fields the current is dominated by edge-like states which 
are not affected by surface disorder. Semiclassically one 
can view carriers motion as short orbits bouncing off the 
same boundary. The associated quantum states have chi- 
rality and are thus commonly refer to as chiral states or 
edge states. The stepwise decrease of the magnetocon- 
ductance observed in regular and chaotic cavities with 
surface disorder is a consequence of the overall depop- 
ulation of Landau levels. It is interesting to note that 
the cavity with bulk disorder shows a markedly differ- 
ent behavior as this type of disorder can, instead, scatter 
carriers between opposite sides of the cavity. The results 
shown in Fig. ^ illustrate the only difference we have 
found up to now between cavities with surface impurities 
or with a number of bulk impurities proportional to L 
JTsfl . Appart from this difference the two behave much 
alike and in line with what one expects to be the behavior 
of quantum chaotic cavities || . 

We have investigated how this steady increase of the 
magnetoconductance is affected by the leads width. In 
order to reduce fluctuations, which are particularly im- 



portant at small W we have averaged the conduc- 
tance over disorder realizations (600 realizations were in- 
cluded in the calculations) . The results for a cavity with 
rather narrow leads attached at the dots in three differ- 
ent ways are illustrated in Fig. ||. At small fluxes (typ- 
ically below 1-2 flux quanta) the expected Lorentzian 
peak characteristic of chaotic cavities || is obtained (not 
clearly visible in the Figure). At higher fluxes, a range 
over which the conductance behaves in a way that is 
strongly dependent on leads configuration, is observed. 
Beyond, the conductance increases linearly with a slope 
which is very similar in the three cases shown in the Fig- 
ure. The crossover to the linear behavior occurs at a flux 
for which the radius of the classical cyclotron orbit r c is 
roughly L/2. In our units (flux and conductance quanta 
$0 = Go = 1 and energy h 2 /(2ma 2 ) = t = 1, where t is 
the hopping integral), r c is given by, 
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where Tiv(E) = (2ysm 2 k x + sm 2 k y )E- At the energy 
chosen here hv(E) w 2.2 Then the flux at which 

r c = L/2 is, for the size of Fig. [| $ = 17<& , which is 
very close to the flux at which the mentioned crossover 
occurs. 

Fig. H shows the averaged magnetoconductance for 
cavities of linear size L = 47 and leads widths in the 
range W=4-20. It is noted that the linear behavior dis- 
cussed above appears in all cases at roughly the same 
magnetic flux indicating that it is only related to the 
cavity size L, as suggested by the discussion above. In- 
stead the slope increases with the leads width. At small 
W (systems with a low conductance) the conductance 
increases very slowly as a function of the magnetic flux. 
For sufficiently large W the conductance reaches a max- 
imum and then decreases stepwise. As remarked above, 
the latter is a consequence of the overall depopulation 
of transversal modes (or Landau levels). The slope of 
the linear part of the magnetoconductance is plotted as 
a function of leads width in Fig. ||. The results can be 
accurately fitted by a W 2 law. The increase of the mag- 
netoconductance can be understood in terms of the in- 
crease of the transmission probability of the transversal 
modes as their edge-like character increases and, conse- 
quently, its sensitivity to surface disorder is reduced. The 
steady increase in the conductance takes place until the 
mentioned depopulation begins to reduce the number of 
modes that participate in the current. Although this ar- 
gument seems plausible, it cannot explain quantitative 
features of the results such as the linear relation between 
the conductance and the flux or the increase of the slope 
as the square of the leads width. This issue is addressed 
in the following Section. 
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IV. DISCUSSION 

The results discussed in the previous subsection 
ressemble those predicted by Beenakker and van Houten 
for the magnetoconductance of two point contact in sc- 
ries fief and the related experiments of Staring et al 
P). Under the hypothesis that transmission between 
point contacts occurs with intervening equilibration of 
the current-carrying edge states, the authors of [[[o| de- 
rived the following expression for the conductance (in 
the following we take the conductance and flux quanta 
Gq = $o = 1 and do not include spin degeneracy), 
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where Ni are the number of occupied subbands in the 
two contacts or leads (i — 1,2) and in the region between 
the contacts, or in the present case in the dot (i = L). 
Disregarding discretness |lC( ], Ni can be written as, 
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where the function /(&) is, 
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with £j = Zj/2r c , U being a characteristic linear dimension 
in the three regions, in the present case, Zj = W\, W^,L. 

We have used Eqs. @) and ( |Io|) to fit the numerical 
results of Fig. ||. We took as fitting parameter the den- 
sity or the Fermi velocity, %v — 2y2wn, and assumed 
the same density in the leads and dot. As shown in 
Fig. [| a satisfactory fitting is obtained for Tiv = 3.65, 
almost twice the actual Fermi velocity in our model (see 
above). The theory reproduces the three regions that 
characterize our numerical results: an almost constant G 
for small flux or r c > L/2, a steadily increasing G up to 
r c « W/2 followed by a steep decrease at higher fluxes. 
It is interesting to note that the theory of Ref. |l(J show a 
better agreement with the numerical results for the case 
in which the leads are attached at opposite corners of 
the dot, than with those for the other two lead config- 
urations of Fig. H A possible reason for this behavior 
relies upon the equilibration assumption in Beenhakker 
and van Houten theory. Equilibration is most likely when 
the leads are not facing each other, as is the case of leads 
attached to opposite corners. Instead, when leads are at- 
tached at contiguous corners direct transmission is more 
probable and equilibration requires higher magnetic fields 
to take place. This is a pictorial illustration of the as- 
sumptions under which the theory of Ref. Jl(| holds. 

In order to check whether a linear relationship between 
the conductance and the flux in a rather wide range of 



fluxes, as indicated by the fittings of Fig.[|, can be un- 
derstood in terms of this theory, we have expanded the 
conductance for small £ and £l > 1- The result for leads 
having the same width and the same Fermi velocity (or 
density) in the leads and dot, is 
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As the maximum in G occurs for £ slightly smaller than 
unity, checking whether G varies linearly with the mag- 
netic field below the maximum requires only to calculate 
the ratio between the coefficients of the second and third 
power of £ in Eq.(|l2|). This ratio is 2.82, indicating that 
the linear term dominates in agreement with our numer- 
ical results. This equation also shows that the slope of 
the straight line is proportional to W 2 . To make a quan- 
titative comparison with the result of Fig. || we rewrite 
Eq. dl2| ) introducing the actual expression for £; the re- 
sult is, 
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The coefficient of the linear term results to be 1.44 x 
lO -4 !^ 2 not too far from the numerical result of Fig. |5[ 
Before ending it is worth to comment on several dif- 
ferences between our model calculation and the theory 
of Ref. 0. We first note that we have assumed that 
the magnetic field is zero outside the dot, which may 
at first sight invalidate the use of Eq. (10). Neverthe- 
less, the Green function of the whole system is calculated 
through Dyson's equation (see Eq. (||). This means that 
the region of the leads close to the dot is distorted by 
the magnetic field, which seems to be enough to validate 
the calculation of the number of occupied subbands by 
means of E q. (|l0| ) . On the other hand whereas in the the- 
ory of Ref. [|10[) each channel contributes with one quanta 
to the conductance, in our case this contribution is ap- 
proximately halved (remember that we work on chaotic 
cavities, see Ref. ||). This seems, however, to be irrele- 
vant as far as the qualitative behavior of the conductance 
is concerned. 



V. CONCLUDING REMARKS 

Summarizing, we have presented a numerical analy- 
sis of the magnetoconductance of quantum chaotic cav- 
ities in a wide range of magnetic fields. For sufficiently 
open cavities the magnetoconductance increases steadily 
reaching a maximum at a flux proportional to the leads 
width. This steady increase of G agrees with the ex- 
perimental observations reported in [Q. Neither regular 
nor bulk disordered cavities behave in this way. Numeri- 
cal results for the average magnetoconductance indicate 
that, for magnetic fluxes larger than that for which the 
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cyclotron radius is approximately L/2 and smaller than 
the flux at which the mentioned maximum is reached, it 
increases linearly with the magnetic flux $ with a slope 
proportional to the square of the leads width. At higher 
fluxes the conductance decreases stepwise. These results 
admit a satisfactory explanation in terms of the theory 
proposed by Beenhakker and van Houten to interpret the 
experimental results for the magnetoconductance of two 
contacts in series. The fact that our results for small mag- 
netic fluxes (r c > L/2) show a better agreement with the 
theory in the case that the two contacts are attached to 
opposite corners of the dot (and, thus, are not facing each 
other), is related to the stronger equilibration of edge- 
states promoted by this lead configuration with respect 
to the other two geometries explored in this work. 
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FIG. 1. Magnetoconductance versus magnetic flux mul- 
tiplied by the inverse of the dot area (both in units of their 
respective quanta), in dots of linear size L and leads of width 
W (connected at opposite corners) with a similar W/L ratio of 
« 0.65. Results for (L, W) — (97,63) chain line, (197,127) bro- 
ken line and (394,254) continuous line, are shown. The results 
correspond to a single realization of disorder and Fermi energy 
E — — 7r/3. Inset: the flux at which the magnetoconductance 
is maximum is plotted as a function of the linear size of the 
dot L. The fitted straight line is, $ mtra /$o = 0.026 + 0.86L. 
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FIG. 2. Magnetoconductance versus magnetic flux (both 
in units of their respective quanta) in 97 x 97 dots with leads 
of width W = 57 connected at opposite corners of the dot. 
The results correspond to dots with: i) no disorder (broken 
line), ii) Anderson impurities with A=6 (a single realization of 
disorder) placed either on all surface sites but those coinciding 
with the leads entrance sites (continuous line) or 2L impurities 
distributed randomly within the dot (thin continuous line). 
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FIG. 4. Magnetoconductance versus magnetic flux in cav- 
ities of linear size L = 47 and leads of width W — 4 (tri- 
angles), 8 (diamonds), 12 (crosses), 16 (circles) and 20 (+), 
connected at opposite corners of the dot. Anderson impurities 
with A = 6 were placed at all surface sites but those coincid- 
ing with the leads entrance sites. The results correspond to 
an energy E — — n/3 and an average over 600 disorder re- 
alizations. Straight lines were fitted for fluxes above 20 and 
below the flux at which the conductance shows a maximum. 
The slopes of the straight lines are plotted as a function of 
W 2 in Fig. |. 
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FIG. 3. Magnetoconductance versus magnetic flux in cav- 
ities of linear size L — 47 and leads of width 8 connected at: 
contiguous corners (circles), opposite corners (crosses), and, 
one corner and the center of the opposite side (triangles). An- 
derson impurities with A = 6 were placed at all surface sites 
but those coinciding with the leads entrance sites. The re- 
sults correspond to an energy E = —n/3 and an average over 
600 disorder realizations. Straight lines were fitted for fluxes 
above 20. The fitted lines are: crosses G = 1.99 + 0.011$, 
circles G = 2.03 + 0.01$, and triangles G = 2.12 + 0.013$. 
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FIG. 5. Slopes of the lines fitted in Fig. ^] versus the 
square of the leads width W. The fitted straight line is 
-0.0019 + 2 x 10" 4 W 2 . 
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